Geometric Control Over the Motion of Magnetic Domain Walls. 
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We propose a method, which enables precise control of magnetic patterns, relying only on the 
fundamental properties of the wire and the choice of the path in the controlled parameter space but 
not on the rate of motion along this path. Possible experimental realizations of this mechanism are 
discussed. In particular, we show that the domain walls in magnetic nanowires can be translated 
by rotation of the magnetic easy axis, or by applying pulses of magnetic field directed transverse to 
the magnetic easy axis. 



Magnetic patterns such as domain walls, bubbles or 
skyrmions, can serve as information carriers in memory 
and logic devices [ij. To efficiently implement such de- 
vices, precisely controlled motion of the magnetic pat- 
terns must be achieved. Several traditional techniques 
have been employed, including magnetic fields 2], field 
gradients Q and electric currents ^j. These techniques 
create a gradient of the energy, which overcomes the pin- 
ning of the magnetic patterns on defects. For exam- 
ple, displacement of a domain wall (DW) in a magnetic 
nanowire with the easy magnetic direction along the wire 
(Fig. [T]) can be induced by applying a magnetic field pulse 
along the wire axis. The resulting displacement of the 
DW depends sensitively on the amplitude and the dura- 
tion of the pulse. 

In this work, we demonstrate that the position of DW 
can be controlled without a gradient of the magnetic en- 
ergy or voltage. Our approach is conceptually similar to 
applications of the Berry'sgeometric phases in topolog- 
ical quantum computing [5|, and to the geometric me- 
chanical control utilized, for example, for precise posi- 
tioning of satellites 0]. The general idea of geomet- 
ric manipulations can be formulated as follows 0, [I]. 
Consider a system exhibiting a continuous symmetry of 
steady states, such as invariance with respect to transla- 
tion along some coordinate x. The system can be con- 
trolled by parameters = {/ii, ^2 ■ ■ ■ l^n\, variations of 
which affect the state of this system but do not break 
the symmetry. When such parameters are varied along a 
closed contour, the final state is generally different from 
the initial one by a shift 5x along the symmetry axis. 
This shift is purely geometrical, namely, it depends only 
on the choice of the contour in the controlled parame- 
ter space but in the adiabatic approximation does not 
depend on the rate of motion along the contour. The 
geometric shift is given by 

5x=(j)A-dfL, (1) 

where ^ is a vector function that depends on the pa- 
rameter values. A well-known example of the geometric 



control in classical mechanics is a falling cat, who is capa- 
ble to rotate its body into an upright position by periodic 
motion of its body parts. 

Geometric control is robust against many errors in 
preparation of the controlling field pulses, which ex- 
plains its success in applications to quantum comput- 
ing and robotics. It is thus worth to explore the pos- 
sibility of similar manipulations with patterns in ferro- 
magnetic materials. We show below that the Landau- 
Lifshitz- Gilbert (LLG) equations describing the evolu- 
tion of classical magnetic systems and the well known 
equations for slowly moving domain walls under exter- 
nal perturbations (3| , lead to a simple realization of the 
geometric control over magnetic patterns such as 180° 
domain walls in ferromagnetic nanowires. 

The geometric control can be implemented, for exam- 
ple, by periodically varying the transverse magnetic field 
or the weak magnetic anisotropy perpendicular to the 
wire. Unlike applied electric currents or magnetic field 
directed along the wire, periodic variation of these pa- 
rameters results in the shift of the DW by a distance 
determined only by the fundamental properties of the 
magnetic system, namely, the shift is robust with respect 
to fluctuations of the driving fields. However, pinning 
on impurities and the intrinsic Peierls-Nabarro barrier 
destroy translational symmetry, and present signifi- 
cant problem to demonstrate the effect. We estimate the 
strength of a magnetic field pulse that can overcome it. 
To demonstrate the effect, we separately consider two dif- 
ferent directions of the magnetic easy axis with respect 
to the axis of the wire. 

DW in a wire with axial easy direction. Consider a 
one dimensional magnetic wire with strong easy axis 
anisotropy along its axial direction x, and a weak easy 
axis anisotropy in the transverse direction forming an an- 
gle 00 with the y-axis. The static magnetic energy is 

E{x) = / dx{J [{dO/dxf -t- sin^ e{d(l)/dxf)\ 

(2) 

+Ko sin^ e + Ksm^ 9 sm^ [cj) - 0o)}, 
where J, and K denote the exchange constant. 
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FIG. 1: The domain wall between two domains with opposing 
magnetization directions parallel to the easy axis along the 



the longitudinal anisotropy and the weak transverse 
anisotropy strengths, respectively, Q — 9{x) is the an- 
gle between the local magnetization vector and the axial 
direction and is the angle between the projection of 
the magnetization on the yz-plane and the y-axis. We 
do not include the magnetic field or the long range dipo- 
lar interactions. The dynamics of 9{x,t) and (f){x,t) are 
given by the LLG equations, which can be written in the 
Lagrange form 



SL d f 5L 



5F 



5X dt \SXj SX 
where X is either 9{x) or 0(x), and 



0, 



L = / dx{E{x) + {Ms/i)4) cos 9}, 
F = (ofAf J27) / dxif)"^ + 02 sin^ 9}. 



(3) 



(4) 



Here, F is the dissipation functional, 7 is the gyromag- 
netic ratio, a is the Gilbert damping parameter, and Mg 
is the magnetization. The DW (FiglT|) is the soliton so- 
lution of the stationary form of Eq. ([3]) , 



6'(a;) = 2tan-^e^^-«-'/'^, (5) 

where A — ^J~JJKq is the width of the DW, and q is its 
average coordinate. 

We are interested in the DW dynamics induced by 
rotation of the weak anisotropy axis, which can be de- 
scribed by variation of 0o in Eq. (I2|) . The simplest imple- 
mentation of this rotation is to physically rotate the wire 
around its axis. We will discuss more practical implemen- 
tations later. We assume that the rotation rate is much 
slower than the characteristic time scale ^ 10~^ sec of the 
magnetic relaxation, allowing a transition from the origi- 
nal variables 9{x,t) and 4){x,t) to variables q{t) and 4'{t) 
that describe the position of the DW and its angle with 
the y-axis [3]. To derive the effective DW Lagrangian 
and the dissipation functional in these new variables, we 
substitute the DW solution ([5]) into ([4]), and integrate 
over X while assuming that q is time-dependent, yielding 



Ldru = {2M,h)^q + 2if Asin2(0 - ^^,{t)), 



Fdw = 



a Ms ( £_ 
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(6) 



The equations describing the DW dynamics are found by 
inserting Eq. ([6]) into Eq. ^ 



2iv:Asin(2[0-0oW]) 



2AIs 2AMsa 
q-\ 



= 0, (7) 



-A<i)/a. 



(8) 



In the adiabatic limit, the last two terms in Eq. ([7]) are 
neghgible, yielding </) ^ (/'o(t). The coordinate of the DW 
is then determined from Eq. ^ 



Q = go - / —dcj^o = go - 27rnA/a, 
J a 



(9) 



where n is the number of rotations of the wire around 
its axis, and go is the initial position of the DW. The 
DW motion thus involves simultaneous rotation together 
with the anisotropy axis, and shift along the wire. Eq. ^ 
directly corresponds to the general geometric control ex- 
pression IT]), where jl is represented by the cyclically var- 
ied Cartesian components of the weak anisotropy axis, 
and Afi,g = — A/a is the connection that relates the evo- 
lution in the control parameter space with the evolution 
in the "fiber space" represented by the DW coordinate q. 
For a steady rotation of the wire with angular frequency 
uj, we find q{t) — qo — i-S- the DW moves with a con- 
stant velocity v — —a; A/a. Equation ^ shows that the 
effect is geometrical in the sense that the displacement 
of the DW depends only on the change in the control 
parameter (jjQ , but does not depend on the rate of change 
as long as the adiabatic approximation is valid. Note 
also that the velocity does not depend on the strength of 
transverse anisotropy K. Such robustness with respect 
to the properties of the wire and the rate of the parame- 
ter variation is a general feature of the geometric control 
making it promising for practical applications. 

It can be shown that the rotation of the wire can over- 
come the DW pinning as long as u > Hc^, where He 
is the critical axial magnetic field required to move the 
DW. In permalloy Py^NigoFeao, i?c ~ 1-5 G 0, yield- 
ing I' ~ uj/{2tt) = 4 MHz, which is too large for practical 
implementations. However, the equations of the DW mo- 
tion derived above remain valid for any other mechanism 
driving the rotation of the DW around the wire axis. We 
note that the perturbation does not need to be steady. 
For example, one can use a short pulse of the circulating 
transverse magnetic field 



Hy = Hoy{t) cos{ujt), = Ho^{t) sin(wi). 



(10) 



with amplitude sufficient to rotate the magnetization in 
the DW, and frequency 1/ > A MHz. For the DW size 
A ^ 20 nm and a ~ 0.01 [ll|, a pulse of duration 
T ^ 1 lis can move the DW by a measurable distance 
27rri/A/a ~ 0.06 mm. The main challenge for the pro- 
posed experiment will be to create an rf field sufficient 
to overcome the transverse anisotropy of the nanowire, 
which for Py is dominated by the shape anisotropy. 
Standard nanofabrication techniques can produce ap- 
proximately square Py nanowires with a cross-section 
50 X 50 nm, and relative variations of the transverse di- 
mensions not exceeding 20%. To overcome the result- 
ing shape-induced anisotropy in Py then requires exper- 
imentally attainable Hqz w 100 Oe ij] and Hoy - pulse 
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FIG. 2: The domain wall between two domains with opposing 
magnetization directions parallel to the easy axis transverse 
to the wire. 



amplitude along transverse anisotropy axis can be much 
smaller. 

DW in a wire with transverse easy direction. Consider 
a ID wire with a strong hard axis along the wire, and 
weak anisotropy in the transverse direction, described 
by the energy functional ^ with i^o < 0, if > 0, and 
A'ol > K. This anisotropy favors the magnetization di- 
rection in the transverse plane, as shown in Fig. [51 If it 
is sufficiently strong {\Kq\ ^ K), then 9 = tt/2 at any 
location along the wire. The energy and the dissipation 
functionals then can be simplified 



dx{J{d(t)/dx)^ + K sm^{(j) - ^o)}, 
F J dx{aMs^^/2-/}. 



(11) 



Without the driving fields (0o ~ const) , the ground state 
is doubly degenerate at = 0o or = 0o + tt. The DW 
solution connecting these two states is given by 



^'^"'(x; 00, q)=(bo + 2 tan"! e(-?)/^. 



(12) 



The motion of this wall under circulating transverse mag- 
netic field was predicted in [l^. We will explore it from 
the point of view of the geometric control Eq. ^ . The 
evolution equation of the angle (f>{x) has a relaxation form 



1 



dt(f> = 2Jdl(j) -Ksin[2{(j) - 0o)]- 



(13) 



Assume that is slowly changed in infinitesimal steps 
(5(/)o, allowing the system to relax after each such step. 
Assume that at some moment of time, the position of 
the DW is g, and the anisotropy direction is rotated 
from 00 ~ to 00. Immediately after this rotation, 
the magnetic state ceases to be the stationary solu- 
tion of Eq. 1131 To find the new equilibrium state to 
which the DW will now relax, we linearize the evolu- 
tion operator comprising the right-hand side of Eq. 
around the equilibrium solution ([12]), taking the form 
b = 2Jd^ldx^ - 2ifcos[2(0'*"'(x;0o,g) - 0o)]. The DW 
can be expanded in terms of the eigenstates Un(x) of this 
operator as Q 

0'*'"(x; 00 - -^00, g) = <f'"{x- 00, q)+Y. GnUnix), (14) 

n 

where C„ are time-dependent coefficients characterizing 
the evolution of the DW. Due to the relaxation, all the 
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FIG. 3: Evolution of the profile of the DW ^ 
parameters used in the calculation are 
uj = -0.00035, K^l,J = l,Hy=0.~' 



x,t) - uit. The 
= 1, 00 = i^t, 
Time t and coordinate 



X are in arbitrary units. 



coefficients C„ will eventually decay to zero, except for 
Co which corresponds to the zeroth mode of the operator 
b, Uo{x) = a,0'^'"(a:;0o,g) = A/ cosh [(a: - q)/A] ^. 
Multiplying ()14p by uo(x), then integrating over a: and 
using the property dxuo{x)(j>'^'^ {x; 0o, 9) = 0, we find 



Co 



^'^'?^o /_!^ dx{uo{x)d^g(j)'^'" {x; 0o, g)} -7rA((50o) 



/-^^2;M2(a;) 



(15) 

The new equilibrium state (j){x) = 0'^"'(x; 0o, q)+CoUo{x) 
is actually a new DW solution shifted by = — Cq. This 
can be shown by expanding the shifted DW solution as 
0'*"'(x; 00, g+5g) « 0''"''(x; 00, g)-a.0'*'"(x; 0o, g)<5(?, and 
noting that uo{x) — dx(j)'''^{x; 0o, <?). For cyclic evolution 
of the control parameter, the geometric shift now reads 



5q 



A^„d(j)o = nn^A, A^^ = 7rA/2. 



(16) 



In contrast to the case of wires with axial easy direction, 
the shift is independent of the Gilbert damping a as long 
as the perturbation is adiabatic. Fig.[3|shows the numer- 
ical solution of Eq. ([T3|) demonstrating the DW motion 
consistent with Eq. ([THj). 

The forces driving the DW due to the rotation are 
suppressed by a small parameter 0ot, where r ~ 10~® — 
10~^s is typical magnetization relaxation time. They 
are considerably smaller in this geometry than pinning 
forces. Hence, rotation alone is not sufficient. However, 
in combination with a strong periodic perturbation suf- 
ficient to overcome the pinning, the geometric shift can 
become important. The strong perturbation can pro- 
duce a zero effect on average, while the small geometric 
shifts accumulate into a significant net displacement. A 
well-known analogy of this configuration in classical me- 
chanics is the Foucault pendulum. Its fast dynamic oscil- 
lations average to zero, while the geometric angle caused 
by the Earth rotation slowly accumulates. 
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FIG. 4: Evolution of the DW profile for the same parameters 
as in Fig. [3] in an additional magnetic field MsHy = 0.012. 



An example of such a perturbation is provided by a 
fixed transverse magnetic field Hy , which is much weaker 
than the effective transverse anisotropy field, but suf- 
ficient to overcome the DW pinning. We consider the 
same rotating wire, whose energy (jlip is now modified 
according to 



E' = E- HyMs / dxcos(0) 



(17) 



Solving Eq. [3] with a constant 
induced DW velocity 



we obtain the field- 



7-ffj,Acos(0o)/a. 



(18) 



In the rotating wire frame, the magnetic field becomes 
oscillating. According to Eq. [THl it moves the DW in 
one direction during one half of the rotation cycle, and 
in the opposite direction during the other half. These 
oscillations are superimposed on a weak but accumulat- 
ing geometric effect due to the rotation of the anisotropy 
axis. Fig. [5] shows the results of simulations for the same 
system as in Fig. |31 but in the presence of the transverse 
magnetic field. DW oscillation due to the magnetic field 
is dominant during each rotation period, but the accumu- 
lation of geometric shifts becomes apparent after several 
full rotations. This effect can possibly contribute to pu- 
rification of magnetic samples from topological defects 
(not only domain walls) during sufficiently long rotation 
in an external magnetic field. 

In summary, we demonstrated that the problem of a 
classical dissipative domain wall motion contains a use- 
ful geometric structure enabling new schemes for robust 
magnetic manipulations. We discussed two examples of 
the DW geometric shift in ferromagnetic wires with dif- 
ferent directions of magnetic anisotropy. We showed that 
adding a periodic perturbation can alleviate the poten- 
tially detrimental effect of pinning on impurities. The 
concept of geometric control is not restricted to the two 



proposed models, and our work should stimulate efforts 
to find its practical implementations in other magnetic 
materials. Other magnetic patterns, such as spin-fiop 
states, spiral waves, and skyrmions, may lend them- 
selves to similar geometric manipulation schemes. Os- 
cillating perturbations have been studied in many sys- 
tems described by the Ginzburg-Landau type of equa- 
tions, includin g d omain walls motion in liquid crystals 
[3, [3, 0, 0, 12, 18 1 . It may be possible to similarly 
interpret many of these results in geometric terms, as 
expressed by Eq. ([1]). Finally, it would be interesting 
to explore the geometric domain wall shifts in systems 
with other symmetries of the order parameter, e.g. in 
ferroelectrics and antiferromagnets. 
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